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Abstract. - The problem of ring formation in solutions of cylindrical micelles is reinvestigated 
theoretically, taking into account a finite bending rigidity of the self-assembled linear objects. 
Transitions between three regimes are found when the scission energy is sufficiently large. At 
very low densities only spherical and very short, rod-like micelles form. Beyond a critical density, 
mainly rings but also rod-like chains appear in (virtually) fixed relative amounts. Above a 
second transition both the length of the linear chains and the relative amount of material taken 
up by them increase rapidly with increasing concentration. The mass accumulated into long, 
semi-flexible worms then overwhelms that in rings. The ring-dominated regime is very narrow 
for semi-flexible chains, confirming that the presence of rings may be difficult to observe in many 
micellar systems, and indeed disappears completely for sufficiently low scission energy and/or 
large persistence length. 



Solutions of highly elongated, cylindrical micelles are arguably among the best studied 
of the so-called equilibrium polymeric systems [Q. Equilibrium polymers are formed in a 
reversible polymerisation process and therefore are in chemical equilibrium with each other 
- monomeric material is continually exchanged between the assemblies. An aspect not at all 
well understood is why ring closure seems to be unimportant in solutions of linear micelles, 
although this would remove unfavourable free ends ("end caps") from the solution. (Closed 
loops have been observed in electron microscopic images of linear micelles Q , but apparently 
they occur in too low concentrations to significantly influence the properties of micellar systems 
[0.) In other equilibrium polymeric systems, such as liquid sulphur, the presence of rings is 
on the other hand thought to be all-important |J. According to mean- field theory ||, rings 
must indeed overwhelmingly dominate the aggregate population in solutions of self-assembled, 
flexible polymers, but only when the monomer fugacity is below some threshold value. At 
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fugacities above this threshold, the concentration of rings should level off to a constant value. 
Linear chains then form an increasingly important fraction of the aggregated material. In the 
limit where the end caps carry an infinite free energy penalty, the crossover becomes sharp and 
is reminiscent of a phase transition, sometimes referred to as the polymerisation transition 

Ignoring the subtle effects of excluded volume interactions which give rise to departure from 
mean-field behaviour j| ||, [|, ^j, the self-assembled linear chains behave as if the rings were 
not present above the polymerisation transition, although the rings do deplete a fixed amount 
of material available to linear aggregation |?J . The amount of material incorporated into rings 
is then roughly equal to the volume fraction at the polymerisation threshold which, for flexible 
aggregates, is predicted to take place at an estimated value of the order ten per cent. This is 
at variance with experimental observations on micellar solutions, where in some cases giant, 
polymer-like assemblies are known to arise at volume fractions as low as one tenth of a per cent 
[0. Some workers have conjectured that a finite bending rigidity of the micelles could suppress 
rings smaller than a certain size, thereby reducing the contribution of rings [Q. Others argued 
that the formation of rings is anyway more strongly attenuated than any estimate based on 
a Gaussian bead-spring model can predict. This is essentially because the number of places 
a continuous ring can break is much greater than that of a discrete ring |t| . Although both 
views are not entirely inaccurate, we shall see below that the issue is decidedly more complex 
and worth of study. 

In this Letter, we reinvestigate linear self-assembly in three spatial dimensions by explicitly 
taking into account the finite bending rigidity of the chains. We find three regimes, one 
where spherical and possibly also short, rod- like aggregates are overwhelmingly present, one 
where closed rings constitute a significant if not dominant portion of the micelle population, 
and one where long, worm-like assemblies emerge to become the dominant species, see figure 
1. For semi-flexible linear assemblies containing many monomers per persistence length, the 
crossovers between the various regimes occur at such low volume fractions that for most 
practical purposes ring formation can be ignored. This explains experimental observation, 
and is in line with an earlier (less elaborate) analysis of Porte 0. We also find conditions 
under which the ring phase disappears altogether. 

We consider a micellar solution above the critical micelle concentration, cmc. We assume 
that at (or slightly above) the cmc the micelles are essentially spherical, and that we can 
disregard the presumably constant amount of the dissolved surfactant not incorporated into 
micelles. In our model, the spherical micelles are viewed as the basic monomeric units that 
self-assemble into elongated, continuous chains. Whether in actual fact single surfactant 
molecules self-assemble into cylinders or whole groups of surfactants is inessential to the theory. 
The cylindrical assemblies can be open (linear) or closed (circular), and are assumed to be 
semi-flexible. This implies that the length each "monomer" contributes to the total length is 
considerably smaller than the persistence length. Self-interactions of the excluded volume type 
are comparatively weak in solutions of semi-flexible chains so the chains may be treated 
as ideal. The dimensionless number density of assemblies of aggregation number N is denoted 
p(N), where the subscripts and 1 shall be used to refer to closed rings and linear chains, 
respectively. Balancing the loss of translational entropy of a monomer when incorporated into 
an aggregate against the free energy gain of association immediately gives, in the mean field, 
p(N) = Z(N) exp(/}iV) for the equilibrium size distribution. The chain partition function 
Z(N > 1) depends on the topology of the aggregate as well as on the prevailing flexibility 
mechanism. We take as a reference the "internal" partition function of a free monomer, 
and define Z(\) = 1. The (dimensionless) chemical potential of the monomers, fj,, can be 
related to the total volume fraction of aggregating material, 0, via the normalisation condition 
4> — 4>q + 4>i, where 4>q = J2n=2 N po(N) and <fii = Y^n=i ^ PiO^O- I n our calculations below, 
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Fig. 1. - Schematic diagram of the aggregation states of semi-flexible micelles. Volume fraction 
4> vs the scission energy E. (i) Monomeric phase (spherical micelles), (ii) ring phase consisting 
mainly of rings and a few rods, both about a persistence length long, (iii) worm phase of long, 
worm-like aggregates, many persistence lengths longs, and a few relatively short rings. Here I 
is a dimensionless persistence length and unimportant prefactors of order unity have been left 
out. 

Fig. 2. - As figure 1, now given as a function of the scission energy E, and the chemical 
potential fj, times the dimensionless persistence length I. Indicated are also the cross-over 
chemical potentials /i* and /i**, discussed in the main text. 



we replace summations by integrations whenever the mean aggregation numbers are large, 
that is, whenever (N) = <fo/E~=a Po(N) > 1 and/or (N) 1 = 0i/£™ =x pi(N) > 1. 

Given that the size distribution is directly proportional to the partition function Z, one 
would next like to obtain accurate expressions for the partition functions Zo and Z\ for 
arbitrary aggregation numbers N. Within the rigid Gaussian chain model of semi-flexible 
linear chains, discussed recently by Marques and Fredrickson M, the latter can be evaluated 
exactly in the entire range from the rod to the worm limit. We only quote the limiting results 
for the rod and worm limits, which can be expressed as Z\(N) ~ z N ~ 1 ((N — l)/l) 3 ^ N ^ 2 ^ 2 
and Zi(N) ~ z N ~ 1 (2iy 3 ^ 2 in these two limits I > N > 2 and N » /. Here I > 1 denotes the 
average number of monomers per persistence length, presumed large, and z = Xil~ 3 ^ 2 exp(E) 
is a function of I and E, with Ai a constant of order unity. A Boltzmann term associated 
with the (dimensionless) free energy cost of chain scission, E > 0, has been explicitly added, 
penalising the free ends of the open chains. For long, polymer-like micelles values of E up to 
about 20 (ksT) are not uncommon Q, although scission energies as high as 38 (fceT) have been 
reported in the literature M]. Note i) that the free energy of a single chain, — In Z, is extensive 
in the number of monomers N in the chain, as it should, but ii) that there are logarithmic 
corrections near the rod limit. These corrections very strongly suppress rod-like aggregates 
considerably smaller than one persistence length yet consisting of more then a few aggregating 
units. (We shall come back to this later.) As Z X (N) < z Ar " 1 (2/)" 3 / 2 for 2 < N < I, we 
simply set Z X (N) = for the entire range 2 < N < /, and Z X (N) = z Ar " 1 (2/)" 3 / 2 for JV > Z . 
The free monomers are from now on understood to also include short, rod-like micelles a few 
micelle diameters long. This merely renormalises the chemical potential. 

Exact results for the partition function of closed loops of arbitrary flexibility are, as far as 
we are aware, not available in the literature, although limiting expressions have been obtained 
within several model descriptions. The problem of ring closure in the tight bending limit was 
addressed by Shimada and Yamakawa at the level of the worm- like chain model O]. Their 
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result can be written as Zq oc exp(— 2ir 2 l/N), valid when / > JV > 1. (Here we ignore 
additional non-exponential dependences on N/l for reasons of brevity.) In other words, tight 
rings are exponentially suppressed. Again, we arbitrarily set the cut-off at N — I, and put 
Zq = for all 1 < N < I. (There is obviously no monomeric ring.) In the opposite limit of 
long chains, N ^> I ^> 1, a simple scaling argument allows us to deduce the partition function 
of the closed configurations, Zq, from that of the linear chains, Z±. The argument is a slight 
modification of the one given by Porte Q (and others || |4|) for the case of flexible chains, and 
hinges on a consideration of the ratio Zi/Z . The ratio Zi/Z is equal to the probability of 
opening a loop, which must be proportional to i) the number of places the ring can break, N, ii) 
the volume that two neighbouring segments can explore after being disconnected, proportional 
to /V 3 / 2 / 3 / 2 [Q, iii) the increased angular phase volume of these new end segments, I 1 / 2 , and 
iv) a Boltzmann weight, exp(_E), for the creation of the two end caps. Using the previously 
given expression for Z\, we find Zq ~ \qz n N~ 5 ^ 2 l~ 2 for N I 3> 1, with Ao an unknown 
constant presumably close to unity. (Pfeuty and co-workers S set Ao = 0.1 in their analysis 
of self-assembled, flexible chains.) We extrapolate this expression in our calculations down to 
N = I. Our estimate for Zq{N > I) is smaller by a factor of l~ x l 2 than the previous estimate 
of Ref. 0, due to the reduced angular phase volume of the two ends upon ring closure, not 
taken into account in that work. 

The size distributions we thus obtain arc 



if l<iV<2 

AoiV" 5 / 2 /" 2 expOuAT) if N > I () 



for the rings, and 




if TV = 1 

px(N) ~ { if 2 < N < I (2) 

if N > I 

for the linear chains, where in both expressions we have absorbed a lnz in the chemical 
potential, fj, = p, + lnz, and a constant into the scission energy, E = E- ln(2 3 / 2 /Ai), which 
we regard as the actual, measurable scission energy. Using eqs (1) and (2), we evaluate the 
relative contributions of chains and rings to the total volume fraction material <f>, and calculate 
their mean lengths (TV) . Careful inspection of the various contributions to 

4> ^ 2A r 5/2 (expQd) - V^^erfc(V^)) (3) 

and to 

0i ~ z- 1 exp( M ) + I 2 expfal - E) (j^y - (4) 

reveals that, depending on the dimensionless chemical potential fi < 0, three regimes can be 
distinguished. See figures 1 and 2. When E > ln(/ 9 / 2 A ^ we find, at a fixed scission energy 
E, transitions between all three regimes upon increasing the chemical potential: 

(i) The monomer-dominated regime: fj, < //*, where /i» ~ — l^ 1 In (Ai/ 1 / 2 + \q\iI~ 4 
x exp(iS)). For volume fractions <f> < A^f 1 / 3 / 2 cxp(— E), almost all of the material resides in 
"monomers" with (TV), = 1, i.e., in spheres and (possibly) in very short rods. 

(ii) The ring-dominated regime: /i* < fi < where /x** ~ — Z 5 / 4 (2Ao) -1 ^ 2 exp(— E/2). 
When < /i < — Z _1 , rings and chains of about a persistence length long occur in a ratio 
almost independent of the concentration, <pi/4> = 1 — {4>o/4>) — 1/(1 + \qI~ 9 ^ 2 exp(E)). The 
mean lengths (N) ~ (N) 1 ~ Z (1 — (/iZ) -1 ) 3> 1 increase only logarithmically with <fi. When 
— l^ 1 < p, < /i** the mass fraction of worms grows with increasing chemical potential. 
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(iii) The worm-dominated regime: /i** < fj, < 0. The fraction of worms overwhelms the 
ring fraction when /i > /it**, that is, when (f> > 4Ao/~ 5 / 2 . The volume fraction of material held 
in rings levels off to </>o ~ 2\ol~ 5 ^ 2 , as does the mean size {N) ~ 3/. We find, furthermore, 
that in this worm-dominated regime {N) ± = I (cf> — <po) 1 ^ 2 exp(E/2) ^> I. This is consistent 
with Porte's result 0, except for the different scaling of 4> a with 1. 

When E < ln(i 9 ' 2 A^~ ), in other words when the micelles are sufficiently stiff for a fixed 
scission energy E, the ring phase disappears and we are left with a single transition from a 
"monomeric" solution of spherical and short, rod-like micelles to a worm-dominated regime. 
The cross-over takes place when fi > /j,* ~ — exp(— nJ), which for all sensible values of 
I gives — iij, = 0(1) independent of the scission energy. The corresponding cross-over volume 
fraction, however, does depend on E, (j> ~ A]" 1 ? 3 / 2 exp(— E) < A^ 1 ^ 3 -C 1. See also the 
"phase" diagrams given in figure 1 (in ip — E space) and in figure 2 (in E — \i space). Note 
that "monomers" do occur in the polymerised phases (ii) and (iii), albeit only in very low 
concentrations. In the electron micrographs of Ref. [Q, spherical micelles can indeed be seen 
in coexistence with long, polymerlike micelles. 

One may infer from our calculations that the so-called sphere-to-rod transition could ac- 
tually be a manifestation of the two transitions discussed above. This is plausible as the 
intermediate (ring) regime can be very narrow indeed for cylindrical assemblies which are 
quite rigid on the scale of a micelle diameter [Q. Even when the assemblies are not all that 
stiff, the ring regime may be difficult to detect as the worm-dominated phase sets in at very 
low volume fractions of order l~ 5 / 2 . Indeed, for fairly flexible linear micelles, such as those 
formed by the surfactant CPyBr in brine, we (conservatively) estimate that I ~ 4, meaning 
that the worm phase must occur at a volume fraction of the order of I0~ 3 above the cmc. 
(Here I was equated to the ratio of the persistence length and the micelle diameter, which for 
this system are accurately known quantities jL3|.) Obviously, we are stretching the theory to 
the limit of its validity in this example, but it does once more confirm that rings should almost 
always be neglected in micellar systems jjj. 

Perhaps the most surprising conclusion of our analysis is that not only tight rings are 
suppressed, but also rod-like assemblies shorter than, say, a persistence length but longer than 
a few micelle diameters. The reason is that close to the rod limit very few bending modes 
remain that could, with the associated configurational entropy, help to compensate for the loss 
of translational entropy upon aggregation. It has to be noted that the precise impact of the 
freezing out of bending modes on the size distribution appears to be non-universal, i.e. it is 
model dependent. Indeed, other workers [l5|, fl6| do not find the relevant logarithmic terms 
in the single chain free energy at all. (In |fl4|, |f6| this is likely to be due to the approximate 
evaluation of Z\(N), and to the highly unusual choice of chain Hamiltonian in The 
question of how strongly our results depend on the choice of the rigid Gaussian chain model 
needs to be addressed. This particular model exhibits "breathing" modes which are not 
present in the standard persistent or worm-like chain model Q, and could therefore exhibit 
spurious behaviour near the rod limit. It is, however, straightforward to verify that terms not 
identical but similar to those found in the rigid Gaussian chain model do arise in a continuous 
description of the (inextensible) worm-like chain model, when taking the limit N — > oo, I — > oo, 
N/l — > flT7fl . (In this limit the problem can be treated at the level of a Gaussian fluctuation 
theory around the straight-rod configuration.) Again we find that the rod-like configurations 
in pi(iV) are strongly suppressed, be it with a different prefactor of N~ N for 1 <C N <C I. It 
is for this reason that we expect the general features of our "phase" diagram to be universal, 
at least in the limit / 3> I . 

Computer simulation studies of self-assembled lattice chains with a finite bending energy 
have recently appeared in the literature |1S, fill E0|. Unfortunately, comparison with our 
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predictions is not straightforward, for the chains interact strongly in the quoted works (and 
in fact exhibit order-disorder transitions) in contrast to the dilute case studied here. Ring 
formation is furthermore explicitly forbidden in JH], ^0|, where simulations were in addition 
restricted to two spatial dimensions only. (Another difficulty might arise from the rather 
low mean chain sizes of order ten investigated in JlJ, [l9|, ^0|.) Tentatively setting aside 
these reservations, some of the trends obtained from our analysis seem to be preserved in the 
simulations. Indeed, in Es) rings were found to play an insignificant role in the ("disordered") 
worm-dominated phase, while in [ fl9| , ppf an increasing chain rigidity strongly reduces the mean 
aggregation number pl| . This happens, in our calculation, close to the (i) — > (iii) transition 

In conclusion, we have have shown, using a simple mean-field analysis, that a finite chain 
rigidity may be the reason why ring formation is unimportant in many (if not all) micellar 
systems. As mean-field theory is thought to have a very broad range of validity in semi-flexible 
systems due to a combination of weak self-interactions and a narrow scaling regime (§], our 
results should be reasonably accurate all the way from the dilute to the concentrated regime, 
provided the system does not cross over to a liquid crystalline phase where a different behaviour 



has to be expected |18 
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